Influence of spin polarization on resistivity of a two-dimensional 
electron gas in Si MOSFET at metallic densities 
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PACS 72.25.Dc - Spin-polarized transport in semiconductors 

PACS 73 . 40 . Qv - Metal-insulator-semiconductor structures 

PACS 71 .30.+h - Metal-insulator transitions and other electronic transition^ s 

Abstract -Positive magnetoresistance (PMR) of a silicon MOSFET in parallel magnetic fields B 
has been measured at high electron densities n S> where Uc is the critical density of the metal- 
insulator transition (MIT). It turns out that the normalized PMR curves, R{B)/R{0), merge 
together when the field is scaled according to B/B^in) where Be is the field in which electrons 
become fully spin polarized. The values of Be have been calculated from the simple equality 
between the Zeeman splitting energy and the Fermi energy taking into account the experimentally 
measured dependence of the spin susceptibility on the electron density. This extends the range of 
validity of the scaling all the way to a deeply metallic regime far away from MIT. The subsequent 
analysis of PMR for low n> demonstrated that the merging of the initial parts of curves can 
bee achieved only with taking into account the temperature dependence of Be. It is also shown 
that the shape of the PMR curves at strong magnetic fields is affected by a crossover from a purely 
two-dimensional (2D) electron transport to a regime where out-of-plane carrier motion becomes 
important (quasi-three-dimensional regime) . 
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Introduction. — Following the experimental obser- 
vation of the the metal-insulator transition (MIT) in two- 
dimensional (2D) systems (for a review, see, for example, 
Ref. [1]), the resistivity of various 2D systems, including 
silicon-based metal-oxide-semiconductor field-effect tran- 
sistors (MOSFETs) and GaAs-bascd hcterostructures, 
was measured in a parallel magnetic field in the vicin- 
ity of MIT PHS] • These studies were largely motivated by 
the still open question concerning the influence of the elec- 
tron spin polarisation on the MIT in 2D (indeed, a parallel 
magnetic field lifts the spin degeneracy without otherwise 
influencing the electron motion). 

It has been observed that on the metallic side of the 
MIT, where n > (for Si MOSFETs, ric w 0.8 -:- 1.0 • 
10^^ cm~^, depending on the degree of disorder), the sam- 
ple resistance R increases with the parallel magnetic field, 
saturating above some characteristic value of the field. Be. 
This value is interpreted as the field corresponding to com- 
plete spin polarisation of the electrons. Possible mccha- 
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nisms of the influence of the spin polarisation on the re- 
sistance were discussed theoretically in Refs. [Tl lTOlfTl] . 

It was also observed [SIE] that once the resistivity 
and magnetic field are rescaled as R{B,n)/R{0,n) and 
B/Bc{n), the PMR curves measured at different densi- 
ties merge into one. This enables one to find the values 
of Be for higher electron densities where the full spin po- 
larisation cannot be experimentally achieved: Be can be 
considered as a variable parameter, to be determined from 
the best merging. Below, we will use this prescription. 

In turn, the dependence of B^ on electron density has 
attracted significant attention. In Ref. [7], it was shown 
that Be is an almost linear function of n vanishing at a 
sample-independent value n = n~,^K, 0.8 • 10"'^^ cm~^. This 
tendency of Be to vanish at a finite electron density raises 
the question about a possible magnetic instability in a 
strongly correlated 2DEG, resulting in a spontaneous spin 
polarisation at n < n^. Then the 2DEG magnetic sus- 
ceptibility Xi which is proportional to the product g*m* 
of the effective Landc g-factor and the rcnormaliscd effec- 
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tive mass m*, must be singular (divergent) at 71 = n^. 
Interestingly, critical density of suggested ferromagnetic 
instability as found in Ref. [7] for a (100) Si-MOSFET 
is close to Uf. ~ 0.8 • 10^^ cm~^, the latter corresponding 
to the MIT. On the other hand, according to Ref. |14| . 
the observed gradual increase of the product g*m* upon 
approaching Uc from the metallic side does not lead to an 
actual divergence of x '^c- This disagreement is dis- 
cussed in Refs. [T5|fT6]. Thus, the dependence Bc{n) near 
ric remains an open question in the problem of the influ- 
ence of electron spin polarisation on MIT in 2D. 

In the vicinity of ric, the correct experimental val- 
ues of Be can be obtained only via very low tempera- 
ture measurements because Be shows a strong tempera- 
ture dependence [9]. Meanwhile, far from the MIT, at 
n ric, the degree of the spin polarisation is almost 
temperature-independent. In this region, the value of B^ 
can be determined from the ratio between Zeeman split- 
ting Aez = g* (IbB (where is the Bohr magneton) and 
the Fermi energy ep = 2n(T:}T? / ra* VgV^) , where and 
are the spin and valley degeneracy factors. As a result, 
the experimental values of Be obtained at low temper- 
atures correctly represent the T ~ Q case. At these high 
densities, the 2DEG physics should be described by a con- 
ventional Fermi liquid picture; it is therefore instructive to 
compare it to the results obtained at lower densities where 
the complications related to the proximity to the MIT be- 
come pronounced. 

In this Letter, we report measurements of PMR at high 
electron densities n ^ Uc] the subsequent analysis allows 
for determination of Bc{n) and for a meaningful compar- 
ison with the low-density case n > Uc. 

Experimental. — We measure the resistivity of a Si 
MOSFET sample in a magnetic field parallel to the sur- 
face and directed along the current. The slot-gate sample 
was similar to those described in our previous publications 
[l7l[T8]. The electron density was controlled by the gate 
voltage Vg: n = 1.24 • 10" cm-2(Vb - 0.34V). The PMR 
curves were measured at T = 0.3 K by the standard lock- 
in technique. 

The high sensitivity of our setup allows for measuring 
a weak PMR signal at high electron densities up to n = 
2.1 • lO^^cm"^. The actual values of gate voltages applied 
(and hence the densities achieved) are restricted by the 
gate breakdown at high Vg- 

The data used in Figs. S] and [5] were obtained in an 
earlier measurement at T = 50 mK, partly reported in 
Ref. [7]. 

Results and discussion. — Fig. [T] shows the PMR 
curves R{B)/R{0) measured at different values Vg and, 
hence, at different n. One can see that the PMR decreases 
rapidly with increasing density (the corresponding values 

strong enhancement of m* with decreasing electron density 
n was also observed in a (111) Si-MOSFET [T2]. Analysis of these 
data shows |13j that the behaviour of m* (n) is similar to that in a 
(100) Si-MOSFET. 
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Fig. 1: (colour online) Normalised magnetoresistance 
R{B)/R{Q) of a Si MOSFET in a parallel magnetic field. 
Curves 1 to 4 correspond to gate voltages Vg = 4V, Vq = TV, 
Vg ~ 12V, and Vg = 17V, respectively. Resistance R{0) in the 
absence of the field is shown in the inset. 

of i?(0) as a function of n are shown in the inset to Fig. 1). 
The data allow us to calculate B^ from the condition that 
the Zeeman splitting and the Fermi energy are equal: 



When calculating the Fermi energy, it is taken into ac- 
count that i/v = 2 for Si(lOO) MOSFET, and that = I 
in a strong magnetic field. The other parameter that en- 
ters Eq. dD) is the renormalised value x oc {g*m*) which 
depends on the density. This latter dependence was re- 
ported in Ref. [14] in the dimensionless scale {g*m*)/2mh 
vs. Tg. Here, 2 is the bare (^-factor for Si, nih = O.IQtoq 
the "bare" effective mass in Si (100), and rg is the ra- 
tio between the Coulomb interaction energy and Fermi 
energy, = mbe^/47reeofi(7rn)^/^ ( with e = 7.7 the di- 
electric constant of the interface between Si and Si02, 
eo is vacuum permittivity). For Si MOSFET, one finds 
rs = 2.65 • (10^2 cm-2/n)i/2_ The calculated values of B^ 
are shown in Table [T] 

Because of the high electron densities, we are far from 
the saturation regime, and the values of B^ can be tested 
only by merging all the PMR curves together by scaling 
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Table 1: Calculated values of Be- The dependence of 
g*m* /2mb on rs is taken from Ref. [14]. 
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Vg, V 
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g*m* /2mh 


-Be, T 
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0.46 


3.91 


1.89 


24.5 


2 
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0.84 


2.89 


1.64 


41.4 


3 


12 


1.45 


2.20 


1.47 


98.8 


4 


17 


2.08 


1.84 


1.38 


151 



them with B/Bc- Fig. [2] shows the result of this scaling. 
One can see that a good merging of all the curves supports 
our calculated values of B^. Finally, in Fig. [3] we plot the 
dependence of Bc{n). We note that at lower densities, the 
values of Bc{n) obtained from Eq. ([T|) are close to those 
determined directly from the saturation of the PMR at 
high fields dlZllIl]. 

We will now discuss the validity of the Fermi liquid-like 
approach which underlies Eq. ([T]) at these lower densi- 
ties. Fig. m shows the PMR dependencies measured for 
the low electron densities n > ric (1-^3- lO^^cm"^) at 
T = 50mK and partly published in Ref. [?]■ In Fig. [SJ we 
plot these PMR curves i?(B)/i?(0) as functions of B/B^, 
where Bc{n) is obtained from Eq. ([l} with appropriately 
renormalised [14] values of {g*m*). The result is shown 
in FiglUin the inset. We see that at the densities below 
n = 1.68 • 10"cm"^ the PMR curves at low fields devi- 
ate from the common scaled curve describing the higher- 
density regime. We thus conclude that at these lower den- 
sities, our calculation oi Bc{n) [see Eq. ([T])] becomes inac- 
curate. Indeed, the data for g*m* [14] used in our Fig. [3] 
were taken at 300 mK, while the MR curves were measured 
at 50 mK. It is known [5] that in the region near MIT the 
value of Be (or equivalently g*m*) shows strong tempera- 
ture dependence. Hence the data of Ref. [14] might not be 
appropriate for the analysis of our R{B) curves at lower 
n, taken at about 50 mK. However, the higher-density 
regime, where Eq. (jlj works perfectly, provides a reliable 
reference point for evaluating B* (n) from the requirement 
that the initial parts of all the R{B / B*) / R{0) coUapsc. 
The values of B* , obtained in this way, are shown in Fig. 
[3] in the inset (dashed line). One sees that the values of 
B* indeed deviate from that of Be (solid line) at smaller 
n. The PMR curve merging is illustrated in Fig. [5] (main 
panel). One may conclude that for the lower concentra- 
tions, n - 1.68 • 10"cm-2, the value of ksT with T = 300 
mK becomes smaller that the energy scale of the critical 
fluctuations, and therefore the values of g*m* obtained in 
Ref. [14] no longer correspond to the true low-temperature 
limiting case (cf. Ref. [H]). Nevertheless, we note that in 
Ref. [14], the data for g*m* were obtained from the analy- 
sis of Shubnikov-dc Haas oscillations, whereas in our Fig. 
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Fig. 2: (colour online) Normalised magnetoresistance rescaled 
as a function of B/Bc- Curves 1 to 4 correspond to gate volt- 
ages 4V, TV, 12V, and 17V. The position of each curve number 
on the plot marks the end of the corresponding curve. 



[3]the values of i?* (suggesting slightly different g*m*) were 
determined using an entirely different method, viz., the 
merging of all the MR curves. The absence of a precise 
agreement between the two is therefore not surprising. 

It is seen from Fig.|4]that down to n = 1.68 • lO^^cm^^, 
MR curves arc saturated above B = Be, which shows that 
the value of B^ indeed corresponds to the complete spin 
polarisation. The behaviour at lower values of n is some- 
what different. We note that at n < 1.53 ■ lO^^cm"^, 
the MR curves at B = B^ exhibit only a bend followed 
by a further increase of resistance with increasing mag- 
netic field, and this behaviour becomes more pronounced 
with decreasing n. This feature is not related to the pe- 
culiar physics underlying the saturation of resistivity at 
B = Be a,s outlined above. Rather, we find that this fur- 
ther increase of resistance reflects a crossover from 2D to 
a quasi-3D mechanism of conductivity, due to the carrier 
motion in the out-of-plane direction which begins to play 
a role in a strong magnetic field, B ^ i?2D- At B = B2D, 
magnetic length A = {h/eBy^^ becomes comparable to 
the thickness of 2D layer Z, viz. A = aZ. Here, a is 
a pre-factor of the order of unity, to be determined phe- 
nomcnologically; we find that a = 0.4 0.5 yields a good 
fit (see the arrows in Fig. |4]). 

The values of Z{n) (as estimated in Ref. [20]) and of 
X{B) are plotted in Fig. [6] (lower and upper curves, re- 
spectively). The dependence of B2D on n is shown in the 
inset. One can see that within the experimental interval 
of magnetic fields (up to 14 T), conductivity of samples 
with higher n always has a purely 2D character. However 
for smaller n, when the corresponding B2D is below 14T, 
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Fig. 3: (colour online) Be as a function of n, Eq. ([T|, for 
high electron densities. Values given in Table [1] (and used in 
Fig. [2]) are shown by bullets. The inset shows a blow-up of the 
low-density region. 

the measured conductivity exhibits a transition from 2D 
to quasi-3D character with increase of B. This explains 
the absence of saturation and further increase of sample 
resistance R after the bend in Fig. 4. Similar MR be- 
haviour had been observed in bulk (3D) heavily doped Si 
with metallic conductivity |21j . 

Conclusion. — When viewed in the context of the 
earlier work on the lower-density Si MOSFETs, our results 
suggest the following generic conclusions: 

(i) We have found that the PMR is present in the high- 
density 2DEG far away from the MIT. In this region, the 
system is essentially a Fermi-liquid and its transport prop- 
erties (including the PMR) are not directly related to the 
occurrence of the MIT at lower densities. This provides a 
further support to the theoretical studies [IO|fTT] which ob- 
tain the PMR based on the Fermi liquid description. We 
speculate that a related magnetic field dependence would 
likely arise in the diffusion coefficient of a 2DEG without 
impurities. 

(ii) We emphasise nevertheless that the PMR observed 
in the high-density regime must have the same underly- 
ing physics as the much stronger effect observed closer 
to the MIT. This is unequivocally illustrated by the fact 
that the same dimensionless scaling (in terms of B / Bc{n)) 
works throughout both regimes, which connect seamlessly. 
The scaling reflects a Fermi liquid-type renormalisation of 
the effective 2DEG parameters (in this case, g* and m*). 
This suggests that the Fermi-liquid phenomenology itself 
is valid all the way down to the vicinity of the MIT, as 
long as the system remains in the ballistic regime. At 
yet lower densities, in the diffusive regime, the non-trivial 
interplay between the interactions and disorder, much dis- 
cussed in the literature (see, e.g., Ref. [12]), strongly af- 
fects the Fermi liquid parameters, resulting in a strong 
increase of the spin susceptibility with decreasing tempcr- 
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Fig. 4: (colour online) Magnetic field dependence of nor- 
malised magnetoresistance in the region of lower carrier den- 
sities (marked in the units of lO^^cm"^). Arrows denote the 
respective values of B2D calculated for a = 0.45 as discussed 
in the text. 
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Fig. 5: (colour online) Normalised magnetoresistance 
R{B)/R{0) scaled as a function of B/Bc (inset) and as a func- 
tion of B/B* (main panel). Here, Be is given by eq. and 
B* is discussed in the text. Values of electron density n are 
in units of lO^^cm^^; for n > 2.21cm^^, the position of each 
value on the plot marks the end of the corresponding curve. 
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Fig. 6: (colour online) The dependence of conducting layer 
thickness Z on the electron density n (lower scale), and the 
scaled magnetic length qA (for a = 0.45) as a function of the 
magnetic field (upper scale). The resultant dependence of Bid 
on n is shown in the inset. 

ature, but this is beyond the scope of the present paper. 

(iii) We find that with an increase of the parallel mag- 
netic field B beyond a certain value B^d at a fixed n, the 
character of conductivity undergoes a crossover from two- 
dimensional to a quasi-three-dimensional one when the 
magnetic length A becomes small in comparison to the 
layer width Z . This is directly observed in the regime of 
smaller densities rt, where the value of B2D lies within the 
experimentally studied range of B up to 14 T. The quasi- 
3D behaviour at high fields is reflected in a continued in- 
crease of sample resistance with B instead of saturation. 
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